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ABSTRACT 

Luminous blue variable (LBV) stars are very massive, luminous, unstable stars 
that suffer frequent eruptions. In the last few years, these stars have been proposed as 
the direct progenitors of some core-collapse supernovae (SNe), particularly Type Iln 
SNe, in conflict with stellar evolution theory. In this paper we investigate various sce- 
narios wherein LBV stars have been suggested as the immediate progenitors of SNe. 
Many of these suggestions stem from the fact that the SNe appear to be expanding in 
a high density medium, which has been interpreted as resulting from a wind with a 
high mass-loss rate. Others arise due to perceived similarities between the SN charac- 
teristics and those of LBVs. Only in the case of SN 2005gl do we find a valid possibility 
for an LBV-like progenitor. Other scenarios encounter various levels of difficulty The 
evidence that points to LBVs as direct core-collapse SNe progenitors is far from con- 
vincing. High mass-loss rates are often deduced by making assumptions regarding the 
wind parameters, which are contradicted by the results themselves. A high density 
need not necessarily imply a high wind mass-loss rate: wind shocks sweeping up the 
surrounding medium may give a high density shell with a low associated wind mass- 
loss rate. High densities may also arise due to wind chimps, or due to a previous LBV 
phase before the SN explodes as a Wolf-Rayet star. Some Type Iln SNe appear to 
signify more a phase in the life of a SN than a class of SNe, and may arise from more 
than one type of progenitor. A Wolf-Rayet phase that lasts for a few thousand years 
or less could be one of the more probable progenitors of Type Ilns, and channels for 
creating short-lived W-R phases are briefly discussed. 

Key words: circumstellar matter; stars: massive; stars: mass-loss; supernovae: indi- 
vidual: 2005gl; stars: winds, outflows; stars: Wolf-Rayet 



1 INTRODUCTION 

Core-collapse supernovae (SNe) arise from stars which have 
a zero-age main-sequence mass > 8Mq, and most likely 
> IIMq. Surprisingly, the evolution of these stars to- 
wards the SN explosion, and their fate, is not well known 
after years of intensive research. Observational programs 
that set out to unearth the progenitors of SNe have so 
far conclusively detected around a dozen, almo st all of 
which appear to be Type IIP SNe (|Smarttl [20091 . and ref- 
erences therein), with red-supergiant (RSG) progenitors. 
The relationship of the remaining types of core-collapse 
SNe to their progenitors is inferred but not determined. 
It had generally been surmised that stars below about 
30 Mq end their lives as RSGs, whereas those above 
30 Mq becom e Wolf-Rayet (W-R) stars before they ex - 
plode as SNe |Falk fc Arnett) Il977l : iPodsiadlowskH 1 1992}) . 
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SN 1987A showed that blue supergiants (BS Gs) could 
also be SN progenitors jSonneborn et all IT987). although 
it may take a binary companion to cause them to explo de 
(|Morris fc Podsiadlowsklll2007l ; [Podsiadlowski et al.ll2007l ). 

Type Iln SNe are relatively rece nt entrants to the SN 
classification scheme l|Schlegel|[l990l ). They are character- 
ized by narrow emission lines (hence the 'n' designation) on 
a broad base, and often (but not always) show very strong 
X-ray and radio emission. These characteristics suggest an 
interaction with a high-density medium, resulting in the 
high X-ray and radio flux. The high dens ity has furthermore 
been attributed to a high m ass-loss rate l|Chugai et al.|[20o3 : 
iGal-Yam fc Leonardll2009l) That , combined with other con- 
sider a tions dKotak fc Vin]jl2006l: iGal- Yam et"al]|2007l : IVinkl 
|2008| ; ISmitbl l2008bl : iTrundle et al.1 |200S| . 120091 ) has led to 
suggestions that the progenitors of Type Iln SNe may be 
Luminous Blue Variables (LB Vs). An LBV, according to 
iHumphrevs fc Davidson! ll 19941 ). is an "evolved, very lumi- 
nous, unstable hot supergiant which suffers irregular erup- 



2 Dwarkadas 



tions" or, in rare cases "giant eruptions like rj Car". The 
high mass-loss rate (> lO _4 M0yr _1 ) during the eruptions 
is mainly what has prompted the suggestion that LBVs are 
Type Iln progenitors. 

Stellar evolution theorists, on the other hand, have 
placed LBVs as a pos t-main sequence pha s e, but not as 



a final pre-SN ph a se JSchaller et al.l 
Langer et al. 19941; IStothers fe Chinl 



1992 



1996 



Langerl 1 19931; 



Maeder et al.l 



2005; iMaeder fc Mevnetll2008l ). In conventional models of 



stellar evolution theory, the LBV phase may follo w a main- 
sequence O-star phase or a H-ri ch WNL phase (|Crowtherl 
2008; iGrafener fc Hamannl [20081 ) . but is always succeeded 
by a H-poor WN phase or a WC phase. Thus stellar mod- 
els require that the LBV star loses its H envelope, be- 
comes a Wolf-Rayet (W-R) star and then explodes as a 
SN. Therefore they were not traditionally considered as 
progenitors of SNe. The hi gh mass-loss rates deduced for 
some SNe such as 1994 W (jchugai et al.ll2004l ) and 2005gl 
l| Gal- Yam fc Leonard! \W0$ ). as well as an explosive event 
that happ ened two years be fore the actual SN explosion 
in 2006jc (|Folev et alj [2007), have prompted several au- 
thors to question whether LBVs could indeed be the im- 
mediate progenitors of SNe. This has led stellar evolution 
experts to explore whether SNe can explode while in the 
LBV stage. Several attempts are underway ([Hirschi et al.l 
l201Ch but are limited by the paucity of knowledge regarding 
LBV and RSG mass-loss rates, and their metallicity depen- 
dence. Although po ssible scenarios have been put forward 
jHirschi et alJl201Ch , it is fair to say that as of the writing of 
this paper, stellar models have so far not succeeded in having 
a star in the LBV stage end its life in a SN explosion. 

This paper takes a closer look at some of the cases where 
LBVs have been proposed as the immediate progenitor of 
SNe. We explore the evidence that has been offered for LBV 
progenitors, and various reasons why other progenitors have 
been excluded. Our goal is to scrutinize various scenarios in 
detail to determine whether the evidence requires an LBV 
as the progenitor, or whether there is room for alternatives. 
Irrespective of whether stellar evolution theory currently al- 
lows for SNe to explode in the LBV stage, we wish to inves- 
tigate whether the observations require that SNe explode in 
the LBV stage. 

The plan of this paper is as follows: In Sj2]we review pro- 
posed scenarios for an LBV star as the progenitor of a core- 
collapse SNe. We investigate aspects of the proposed sce- 
nario not considered by the authors, and examine whether 
the data are compatible with an LBV progenitor. sum- 
marizes our results and outlines our suggestions, especially 
where Type Iln SNe are concerned. 



2 PROPOSED CASES FOR LBV 
PROGENITORS OF SNE 

2.1 SN 2005gl 

This luminous type Iln SN is located in the nearby galaxy 
NGC 266. Pre-explosion Hubble Space Telescope (HST) im- 
ages of the SN location led lGal-Yam et all j2007) to suggest 
that a very bright point source with luminosity L > 10 6 Lq 
was the possible prog enitor star of this SN. This identifica- 
tion was confirmed bv lGal-Yam fc Leonard! (|2009l . hereafter 



GL09) when they observed that the putative progenitor had 
faded away and was no longer visible in HST images. 

Spectra taken 8 days after discovery (between 8 and 
33 days after explosion) show a narrow component to the 
Ha emission with a velocity of 420 km s _1 , and an inter- 
mediate component with a velocity of about 1500 km s _1 
jGal-Yam fc Leonard! 120091 ) . No broad components repre- 
senting typical SN ejecta velocities in the range of 10 4 km 
s _1 is seen. Spectra on day 58 after discovery however do 
show a broad component with a 10 4 km s _1 velocity. The 
intermediate component is missing. The narrow component 
on day 8 is interpreted by GL09 as the velocity of the un- 
shocked progenitor wind. The velocity is deemed too high 
for red supergiants, and too low for Wolf-Rayet stars, but 
about right for an LBV star. The dense wind slows down 
the SN shock to about 1500 km s~\ whose velocity is re- 
flected in the intermediate component. By day 58, according 
to GL09, the shock has crossed the dense wind and entered a 
region of less dense wind, and the emission is dominated by 
fast, unshocked ejecta. Assuming a steady LBV wind, GL09 
determined a mass-loss rate of about 0.03 Mq yr~\ which 
GL09 contend can only arise from an LBV star. 

The large velocities observed in the day 58 spectrum at- 
test to the fact that the ejecta have sufficient kinetic energy 
to expand with initial velocities > 10 4 km s~\ Difficulties 
with this scenario could be associated with slowing down the 
fast-moving ejecta interacting with the ambient medium to 
about 1500 km s _1 in as little as 8 days, or as much as 33 
days depending on the actual explosion date. We investigate 
various aspects of this scenario below. 

Since the ejecta were expanding freely again on day 58 
in this scenario, we estimate that they would have crossed 
the dense LBV wind in 50 -75 days or less (depending on 
explosion date). By day 8 they were already at 1500 km s _1 , 
and would have slowed down even more, we may assume that 
the ejecta had an average velocity about 1250 km s _1 for 
the 50-75 days (depending on explosion date) that it took 
to cross the LBV wind medium. Since the velocity of the 
LBV wind was about 1/3 that of the SN, it means the wind 
in this phase would have lasted for about 3 times more, or 
about 150-225 days. Given the mass-loss rate of 0.03 Mq 
yr _1 calculated by GL09, this means that the total mass 
ejected in this wind was about 0.0125-0.0185 Mq. This is an 
upper limit, as the wind phase could have conceivably lasted 
less than 50-75 days, and the average shock velocity could 
have been even lower (note the argument in the previous 
paragraph), given that it had slowed down from > 10,000 
to 1500 km s -1 in just 8-33 days. 

If we make the extreme assumption that all the CS mass 
is contained in a thin, dense shell, the properties of the inter- 
action have been shown to depend on t he ratio of the mass 
of th e shell to the mas s of the ejecta (Tenorio-Ta gle et al.l 
Il990l ; ID warkadasl 12005). A ratio of 1 or higher is required 
for the shell to have a significant impact, while a ratio << 1 
would ha ve no impact as the ejecta would just plow through 
the shell. iDwarkadas! (|l997l ) also showed that the distribu- 
tion of the mass, whether it is in a thin shell or spread out, 
is not as important as the amount of material. Therefore, 
in order to decelerate the ejecta completely, as seems likely, 
the ejecta mass involved in the interaction must be smaller 
than or comparable to 0.02 Mq. 

The important question therefore is how much ejecta is 
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interacting with this amount of CSM material. This is dif- 
ficult to say unless we can specify the ejecta mass and its 
distribution, and how the reverse shock expands into this 
ejecta, all of which are unknown. However, we can still make 
some estimates as to whether the scenario proposed is consis- 
tent with the observations. Following the arguments above, 
almost all the kinetic energy of the fraction of the ejecta in- 
teracting with the CSM is lost (the shock velocity decreases 
to 1/7 the original, therefore the total energy decreases to 
1/49 the original). The total kinetic energy of the amount of 
ejecta interacting with the CSM was, at most, 0.5 M ej v 3 2 , 
where v si is the initial velocity of this material, which is 
at least 10,000 km s _1 . The maximum luminosity that can 
be extracted from the radiative shock is 0.5 MVg/v w . Given 
the mass-loss rate M = O.O3Af0yr _1 and the wind velocity 
v w = 420kms _1 deduced by GL09, the maximum luminosity 
that can be extracted immediately (when the shock velocity 
is > 10 4 km s -1 ), is 4.5 xlO 43 ergs s _1 . However, it is un- 
likely that the luminosity will be so high when the shock is 
moving so fast. On the other hand, the maximu m luminosity 
on day 8 (33) is 1.5 xlO 41 ergs s~\ (Note that lSmith et all 
|2010l ) find the wind velocity 34 days after discovery to be 
closer to 600 km s~\ which would reduce these numbers by 
about 2/3) . Since as the shock decelerates the amount of en- 
ergy released would be expected to increase, we may assume 
that the average luminosity for the 8 (33) days was about 
1. xlO 42 ergs s _1 . Over 8 days, the total energy released 
is 6.9 xlO 47 ergs, while over 33 days it would be at most 
2.85 xlO 48 ergs. The latter is probably an overestimate, as 
the shock deceleration would be slower over a larger period. 
Equating that to the total energy available for release gives 
the ejecta mass interacting with the CSM to be about 1-3 
x 1O~ 3 M0, and possibly up to a factor of 4 smaller. If the ini- 
tial velocity is higher this mass increases somewhat, whereas 
if the wind velocity is closer to 600 km s _1 this decreases 
some. 

The ejecta mass is an order of magnitude smaller than 
the mass in the LBV wind it is interacting with. Therefore 
deceleration of the ejecta is not an issue. We next check 
whether this amount of mass in the unshocked CSM is suf- 
ficient to provide the requisite Ha luminosity in the narrow 
line seen by GL09. We assume that the average velocity of 
the shock was about 2500 km s" 1 over the first 8 days after 
discovery, or slightly lower if 33 days after explosion. There- 
fore, the maximum radius of the shock would be between 
about 1.7-7.1 xlO 14 cm. This shock would therefore sweep 
up a mass of the CSM equal to at most 1.6 xlO _2 M0. The 
mass of unshocked CSM remaining is about 2.5 xlO _3 M0 
or larger if the smaller time period is considered. Since the 
wind expanded for at most 225 days at 420 km s _1 , it could 
have reached about 8.1el4 cm (somewhat larger if the ve- 
locity was closer to 600 km s" 1 ). The volume of unshocked 
wind is about 7.3 xlO 44 cm 3 . The luminosity of the Ha 
emission per unit volume can be approximately written as 
Lh q = an^hfc, where a = 1 — 3 x 10~ 13 is the recombina- 
tion coefficient. The observed luminosity of 2.8 xlO 39 ergs 
s _1 in the narrow line gives the required H number density 
of the unshocked wind at that radius to be about 2 x 10 9 
cm~ 3 . We note that this is consistent with the actual num- 
ber density of the unshocked wind just outside the shock. 
Furthermore the total mass of H required is smaller than 
the mass of unshocked wind. Given the many assumptions 



that have gone into this calculation, it appears reasonable 
that the unshocked wind has sufficient density to provide 
the required Ha luminosity. Note that if we take the expan- 
sion time as only 8 days then the density constraint is even 
easier to satisfy. 

One final check to be made is whether the fast moving 
shock would become radiative immediately. We consider a 
shock of 11,000 km s _1 on day 1. The cooling time for such 
a shock would be around 1 day at a density of 10 10 cm -3 , 
which is expected for the assumed numbers very close in 
to the star. Thus it seems likely that the shock would be 
a radiative shock, would be able to radiate the observed 
amount of energy in the appropriate time period, and that 
the ejecta mass would be small enough to be decelerated by 
the small amount of LBV material. 

We note that the above calculations for the maximum 
luminosity were carried out assuming the existence of a 
steady wind, to be consistent with the assumptions of GL09, 
but without any supporting evidence for this assumption. As 
we show later in this paper, in many cases a steady wind is 
not supported by the observations. However, the above re- 
sults could be derived in a different way. A lower limit to 
the available luminosity is set by the observed spectrum on 
day 8, especially the Ha line luminosity. If we assume that 
the Ha luminosity is about a tenth of the total available lu- 
minosity, the other numbers follow accordingly. Therefore, 
even if one removes the steady wind constraint, the com- 
puted ejecta mass would not change much, and would still be 
lower than the calculated LBV mass. Removing the steady 
wind constraint would mean that we would not be able to 
easily calibrate the mass-loss rate. However, given the high 
density, which is required for the shock to be radiative (oth- 
erwise a broad Ha line would be seen on day 8), the remain- 
ing numbers would still be consistent. Also, in order for the 
shock wave to become radiative, a high density must exist 
from day 1 immediately outside the stellar envelope. 

Although the calculations are necessarily approximate, 
and could be off by factors of 2-3, it appears that the inter- 
pretation of an extremely dense medium immediately sur- 
rounding the star may be consistent with most of the ob- 
served facts. The fact that the shock was no longer radiative 
in the spectrum taken on day 58, given that the velocity 
must have been quite close to the initial velocity, suggests 
that the density in the wind into which the ejecta are now 
expanding must be lower by more than 2 orders of magni- 
tude, to make the cooling time at that epoch larger than the 
flow time. If this is a wind from an LBV in quiescence, this 
seems reasonable, as the mass-loss rate will be at least two 
orders of magnitude lower. If this is from a prior phase, say 
a RSG wind, then since its velocity must be a factor of 20 or 
so smaller, this means that the mass-loss rate must be more 
than 3 orders of magnitude smaller, or of order lO _5 M0yr _1 
or lower. 

These calculations corroborate the scenario of the SN 
expanding in a very high density medium surrounding the 
star, with a mass-loss rate of about 0.03 Moyr -1 (assuming 
a steady wind) and a wind velocity between 420 and 600 
km s _1 . GL09 suggest that these parameters are typical of 
LBVs and not found in other stars. While this is true, it 
must be noted that even in LBVs, such high mass-loss rates 
are only postulated for LBVs when they are undergoing an 
eruption, not when they are in a steady state. The detection 
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of a very bright source > 1O 6 L0 in roughly the V-band is 
used by GL09 to support their assertion that this was no 
ordinary star. Does this definitely make it an LBV? Our 
knowledge of stellar evolution and the placement of stars on 
the H-R diagram is incomplete, and it is unclear if every star 
with luminosity > 10 6 Lq need necessarily be classified as an 
LBV star. LBVs are marked by their variability, which has 
not, and probably cannot, be demonstrated in this candidate 
star. Of concern also is the described evolutionary scenario 
that one is forced to accept. In this scenario, a star tran- 
sitions into an LBV-like state, undergoes mass-loss with a 
mass-loss rate seen only in r/ Car type eruptions for about 6 
months or so, and then undergoes core-collapse. Did the star 
transition from a quiescent LBV into an eruptive state, and 
explode during, or immediately after, that state? Did the 
star transition from a different phase (say a red supergiant) 
to the LBV-like state, making it an LBV for only 6 months, 
during which time it underwent an eruption and then ex- 
ploded as a SN? If the explosion of a star as an LBV is by 
itself a problem, the explosion of a star immediately after 
erupting as an LBV, with an extremely high mass-loss rate, 
is an even greater problem for stellar theorists to reconcile. 

Therefore, many questions still remain unanswered. 
Given the available information, and with the above caveats, 
we currently rank SN2005gl as making an acceptable case for 
an LBV-like progenitor. Any alternative must have a high 
density wind (for the given velocity of around 500 km s _1 ) 
and a luminosity > 10 6 Lq in the V band, and satisfy all the 
constraints computed above. 

2.2 Quasi-Periodic Modulations of Radio 
Supernovae 

I Ryder et al.l (|2004T l found modulations in the radio light 
curve of the Type lib SN 2001ig, with a periodicity of 150 
days (altho ugh this is based on only two peaks, with a prob- 
able third). ISoderberg et al.l (|2006T l report similar modula- 
tions in the radio lightcurve of SN 2003bg, at 120 and 300 
days, with a probably third at 600 days. In this case it is 
not clear that they are periodic. ISoderberg et al.l (|2006l ) at- 
tribute these to varia tions in the density of the CS medium. 
iKotak fc Vinkl { 2006) suggest that these density variations 
may be due to LBV stars undergoing S Doradus type varia- 
tions, with the density enhancements specifically due to the 
behavior at the bistability limit, which would lead to density 

jumps of a factor of 4. 

The reason for IKotak fc Vinkl |2006l ) suggesting LBV 
progenitors was to suggest a single star scenario rather 
than more contrived binary star scenarios. According to 
IKotak fc Vinkl | 2006), the LBV scenario was in line with 
both the timescales as well a s the amplitudes . How ever 
modelling of the light curves by ISoderberg et all (|2006T ) has 
shown that density enhancements of factors of 1.8, 1.4 and 
1.2 are required to fit the "bumps" in the light-curve. These 
are relatively small enhancements, and it is not clear why an 
LBV progenito r would need to be in voked to explain them, 
as was done bv lKotak fc Vinkl l| 20061 ). In fact it is somewhat 
of a concern as to whether LBV outbursts would give rise to 
such low-level density enhancements; on the contrary, they 
will presumably lead to much larger changes in density than 
the 20-40% required t o explain some of the lightcurve mod- 
ulations. According to lKotak fc Vinkl |2006T l themselves, the 



bistability jumps lead to a density enhancement of a factor 
of about 4. Mass-lo ss rates and veloci ties derived around the 
bistability jump in IVink et al.l (1 19991 ) suggest perhaps even 
larger, not smaller, density enhancements, and no models 
suggest 10-20% enhancements. It is therefore unlikely that 
the density jumps due to S Doradus variations would match 
the modelled density enhancements. An exception would be 
if both the mass-loss rate and velocity enhancements were 
in the same direction, thus leading to a small enhancement 
in the density, which depen ds on their ratio. It sho uld also 
be noted that the models by ISoderberg et all (|2006l ) are ap- 
proximate, and not based on hydrodynamical calculations. 
As they note, the self-similar solutions they have used are 
not truly applicable in a case where the density is increasing, 
thus leading to some uncertainty in the calculated results. 

Even if bistability jumps were adequate, it is not clear 
that they should point only to an LBV. Although LBVs are 
the only objects in which bistability jumps have actually 
been observed, other massive stars, besides LBVs, coul d also 
experience these bistability jumps. IVink et all j 19991 ) have 
shown that the increase in mass-loss rate at the bistability 
jump is mainly due to radiative acceleration by Fe III, and 
therefore should occ ur in hot s t ar atm ospheres where this 
ion is present. And IVink et al l (2000) have derived mass- 
loss rates for all massive O and B stars assuming that they 
are affected by the bistability jumps. 

An added consideration is that if the enhancements are 
due to changes in the wind parameters, then the wind-wind 
interaction will likely result in the sweeping up of the outer 
wind into a dense shell. This will increase the magnitude of 
the density enhancement, suggesting that the actual change 
in wind parameters may be even smaller, and making LBV 
progenitors even less likely. 

The timescale calculation is uncertain as it relies 
on uncertain velocity determ inations, as pointed out by 
IChevalier fc Soderberd (fcOld ). Furthermore, 2003bg was ini- 
tially classified as a Ic due to lack of hydrogen in the 
spectra. It is not known why a H-rich LBV progenitor 
would give rise to a Ic spectrum. Perhaps the star transi- 
tioned into a W-R star just before explosion (Sj3] explores 
such short W-R phase s ). Ac cording to the calculations of 
IChevalier fc Soderberd (|2010) it would have had to happen 
on order of a quarter century before explosion. If this were 
true, the fast W-R wind would be expected to sweep up the 
preceding LBV material into a thin dense shell, which is not 
seen. Also, if this were true, an LBV would not be the direct 
pro genitor of the SN. 

IChevalier fc Soderberd (|2010l ) suggest that these SNe 
lie in a category called cllb, or compact Type lib, with 
radii of a few times 10 11 cm. They find that the progeni- 
tors of these SNe are less massive and more compact than 
LBVs, further discounting LBV s, and perhaps more su gges- 
tive of Wolf-Rayet (W-R) stars. ISoderberg et all (|2006l ) sug- 
gested a single W-R star progenitor with an average mass- 
loss rate of 3 xlO _4 M0yr _1 . However, this mass-loss rate 
is much higher than th e known mass-l oss rate of any sin- 
gle galactic W-R star l|Crowtheijl2008l ). The combination 
of a high mass-loss rate and a short W-R progenitor phase 
may suggest (see also £|3 . 2 1 b elow) that the SN e xplosion 
occurred in a bina ry system (|Eldridge et al.ll2008l ). In fact 
iRvder et alj d2006|) claim to hav e foun d a companion star, 
although Chevalier fc Soderberd (|2010h argue that the star 
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would have a we ak wind and would not create the pinwheel 
like nebulae that I Ryder et al.l l|2006l ) suggest gives the den- 
sity variations. While certainly true, this does not preclude 
the dynamical effects of the binary companion. It is interest- 
ing to speculate whether the stellar orbit of the companion 
(rather than its wind), especially if it was quite eccentric, 
could affect the mass-loss from the star, leading to both an 
enhanced mass-loss rate and a variation in the mass-loss as 
the orbital radius changes. 

In summary, although we can speculate on the progen- 
itor, we are unable to pinpoint it. While strong arguments 
can be made a gainst S Doradus type variati ons based on 
the modeling by Chevalier fc Soderberd |201(J ). there is also 
some uncertainty in the models themselves. These issues re- 
quire further study and testing via both observations and 
modelling. Although we lean towards this scenario not pro- 
viding sufficient evidence to suggest LBV supernova pro- 
genitors, a final assessment requires more observations than 
are currently available, accompanied by more detailed mod- 
elling. 



2.3 Type Iln SNe 

Many Type Iln SNe show signs of strong CS interaction, in- 
dicating that the SN shock may be encountering a medium of 
high density (ISalamanca et al.l T1998I: [Pastorello et al.ll2002l; 
Salamanca! 120031; IChugai fc Danzigei 20031 ; Chugai et al.l 
20041 ; iBauer et al.l 120081 ) , which has been interpreted as a 
sign of a high wind mass-loss rate. In order to investigate 
this assumption, it is useful to understand how the mass-loss 
rates have been calculated in several cases. If the medium 
into which the SN is expanding is a wind with a constant 



mass- loss rate M and wind velocity v w , then the wind den- 
sity decreases as r~ 2 . The mass swept up by the expanding 
SN shock up to a radius R is the mass of the wind upto that 
radius 



M 3W = MR/v 

w 

The total energy in the swept-up material is 



E = 0.5M sw vi 



(1) 



(2) 



where v s is the shock velocity. Therefore, the kinetic energy 
dissipated per unit time, which is the maximum energy that 
can be extracted from the material over time, or the maxi- 
mum attainable luminosity, is: 



dE/dt - 



v w 



(3) 



shock velocity is constant with tim^ If this is not the case, 
then equation [T] is not valid, and neither are the succeeding 
equations, which of course means that equation [3] cannot 
be used. If it is the answer cannot be relied upon. In gen- 
eral, if the density profile is not at least as steep as r~ 2 , 
application of these equations at a given radius will yield 
a misleadingly high mass-loss rate. Unfortunately, in many 
of the cases that they have been used, as illustrated below, 
the values derived contradict the assumptions of constant 
mass-loss rate and velocity, and an r~ 2 density profile. 

For SN 2006gy. ISmith et al.1 (|201Ch find, by using equa- 
tion |31 that the value of both M and Vcsm is varying in 
time (their Figure 25). This is contrary to the assumptions 
inherent in the equation. Furthermore the SN blast wave 
is assumed to travel within this medium with a constant 
velocity, which is incompatible with the derived result of 
a medium with a varying mass-loss rate and velocity with 
time. It is almost impossible to construct a model in which 
the density parameter w = M/v w varies non-uniformly by 
over a factor of 100 (as the authors find) while still keep- 
ing the velocity of the blast wave roughly constant. Their 
argument hinges on the untrustworthy computation of the 
wind density parameter from equation |3j. As they them- 
selves point out, their assumption that the entire luminosity 
comes from energy radiated behind a r adiative shock is c on- 
trary to their own earlier assumption l|Smith et alj|2007l ) of 
a shell-shocked diffusion model. If a substantial part of the 
luminosity arises from diffusion of energy from a shocked, 
optically thick shell, then this would lower their estimates 
considerably. Because of the inconsistencies in deriving these 
numbers, the total mass of swept-up circumstellar matter 
(CSM) is also in doubt, and the high mass- loss rate ques- 
tionable; 

iMiller et al.l (|201dh use equation[3]to compute the mass- 
loss rate in SN 2008iy. They assume a fixed wind velocity 
(although it is measured only at one epoch) and a constant 
blast wave velocity. They find that the mass-loss rate in- 
creases with time, then decreases. Given the increase in lu- 
minosity over time, they assert that the parameter w must 
have increased over 400 days. However, if the parameter w 
increases with time, the density profile no longer goes as r~ 2 
over the 400 day rise period, the shock velocity should cer- 
tainly not be constant, equation [T] is not valid, and neither 
are the others. Furthermore, although they assume a model 
of clumps in a rarefied wind to explain the late-time emis- 
sion, they do not take into account any luminosity arising 
from shocks within the clumps in their mass-loss rate calcu- 
lations, which could lower their estimates significantly. Thus 
the mass-loss rate is inconsistent with the assumptions used 
to obtain it, and therefore unreliable. 



If the luminosity, shock velocity and the wind velocity are 
known, then this equation can be used to compute the mass- 
loss rate, with perhaps an unknown conversion efficiency 
factor added. Some form of this equation has been used to 
compute the mass-loss rate for several Type Iln SNe, in- 



cluding SN 1997ab (ISalamanca et alJll998l), SN 2005gl, SN 
2006gy (ISmith et alj|2010h. SN 2" o08iy ^Miller et al.1 12010? ) . 
and SN 2006tf jSmith et aljlioosh . The important point to 
note from the above derivation is that this equation is valid 
if, and only if, the wind mass-loss rate and velocity are con- 
stant with time, such that its density goes as r~ 2 , and the 



1 We note here for the sake of completeness that, if the shock is 
radiative, the formula 

L = 2nr 2 pv 3 

can be used to compute the pre-shock density, assuming that all 
the energy is radiated away. This does not require a r~ 2 medium, 
but it does require a radiative shock. However substituting p = 

t-^t~ — to get the same equation as above does require assuming 
a constant mass-loss rate and wind velocity. All of the papers 
discussed here have assumed the medium goes as r — 2 and used 
equation [3] to derive the mass- loss rate 
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iMiller et all (|201CJ ) also compute the mass-loss rate 
by two other method s: (1) from the X-ray luminosity 
ijlmmler fc Kunt j l2005h and (2) from the H a luminosity. 
Unfortunately, the equation used to compute the mass-loss 
rate from the X-ray luminosity makes the same assumption, 
of the density of the medium decreasing as r -2 , which is 
contravened by their results. Furthermore the authors use 
a low X-ray temperature of 10 7 K to compute the cooling 
function at an age of 1.5 years, which is highly unlikely for 
a SN at such a young age, and incompatible with the high 
blast wave velocity of 5000 km s _1 that they have assumed, 
which would give a post-shock temperature of few times 10 s 
K. It is highly unlikely that the spectrum of such a com- 
plex dynamical region could be adequately described by one 
temperature, in any case. They have not taken into account 
the possibility that the X-ray emission could be arising from 
both the forward and reverse shocks. No attempt is made 
to take the SN deceleration into account when calculating 
the radius. Finally, we note that the numerical coefficient 
us ed in this equation is incorrect due to a missing factor of 
4 ijlmmler et alJl2006t ). 

In their computation of the mass- lo ss rate from t he Ha 
luminosity, it is not clear what velocity [Miller et al.l (|201Ch 
refer to when using Vsn- Early in section 4.1 they refer to it 
as the velocity of the blast wave overrunning the CSM. This 
would then make this the velocity of the forward shock. In a 
medium as dense as the one they compute, the s hock wave 
would sweep up the material and decelerate. The IChevalierl 
(|l982h self-similar solution gives the shock velocity as going 
as where the deceleration parameter a < 1. The radius 
would then be given by Rsn = VsNtsN / ol; the authors have 
neglected the parameter a. For the method of computing the 
mass-loss rate from the Halpha luminosity, they claim that 
ri is the "inner radius" corresponding to the position of the 
blast wave, but use the same formula n = VsNtsN, which 
is incompatible with their earlier assumption that this is the 
velocity of the blast wave overrunning the CSM. It was also 
assumed that a radius r-z was the outer radius related to 
the fast-moving forward shock, and that ri >> ri. This is 
puzzling, because the outer shock is the one overrunning the 
CSM, which they earlier referred to as the blast wave. The 
assumption that ri » n at day 711 would be incorrect in 
any model. The IChevalierl l| 19821 ) solution predicts that the 
ratio of the outer to inner radii would be more like 20% in the 
self-similar case. In all cases they have neglected the value 
of a in their computations, which would make the radius at 
each epoch slightly larger. 

In their sce nario to explain the origin of the 400d rise 
time for this SN. lMiller et al.l (|2010T l favor the clumpy wind 
scenario. As mentioned, it is then surprising that they do not 
take into account the luminosity from the radiative shocks 
in the clumps, which wou ld reduce their mass-loss rate esti- 
mates. IMiller et all l|2010h also argue against a wind-bubble 
scenario on the grounds that after peak the thin shell is 
overtaken by the ejecta, and the luminosity is powered by 
interaction with the ambient wind. This is not necessarily 
true. The shell could be thick, and the peak reached while 
the shock is still interacting with the shell. Even when the 
forward shock exits the shell and is interacting with the am- 
bient wind, the reverse shock is still interacting with the shell 
for some time and could power the luminosity , and an ex- 
tremely dense wind may not be needed. Finally. [Miller et al.l 



|2010l ) say that addition of a wind-bubble component re- 
sults in an unnecessary transition from an LBV to a W-R 
star. However, this is exactly the kind of scenario we envis- 
age for many Ilns in i]3.4l Without detailed calculations, a 
wind-bubble scenario cannot be excluded. Due to the vari- 
ous questionable assumptions, their mass-loss estimates are 
untrustworthy, and the LBV scenario remains unproven. 

In the case of SN 2006tf, the mass-loss rate derived by 
ISmith et all ([2008), again using equation [3l exceeds 1 Mq 
yr _1 , an astonishingly high rate. The mass-loss rates and 
wind density parameter are found to decrease gradually with 
time over 10 years, and decrease by a factor of 10 after that. 
Note that a mass-loss rate decreasing outwards, combined 
with a constant wind velocity, gives a density profile that de- 
creases faster than r -2 . This again means that the results are 
inconsistent with the assumptions (constant mass-los s rate 
and velocity) used to derive them. ISmith et al.l |2008) have 
assumed a shell-shocked diffusion model to explain the high 
luminosity, so it seems strange that they would calculate 
the mass-loss rate assuming that the entire luminosity arises 
from a radiative shock with 100% energy conversion (which 
is inconsistent with the constant shock speed assumed). 

The mass-loss rate about 20 years prior to explosion 
is found to decrease by a factor of 10, by analysis of the 
narrow lines. Assuming that the narrow lines arise from the 
CSM leads to the high estimate for the mass-loss rate and 
total mass. A plausible alternative may be that the narrow 
line emission arises from a very small section of the CSM, 
such as a denser clumpy component. The ratio of the narrow 
to broad line velocity of a factor of 10.5 suggests that the 
density contrast between the clumps and the surrounding 
medium would be about a factor of 110. The narrow line 
velocity would not be the velocity of the CSM, but that 
of the shock within the clumps. Then the density of the 
medium is a factor of 100 lower and does not contribute to 
the emission, and the large mass-loss rate, and consequently 
large mass of the progenitor, are not necessary. The narrow 
and broad line components arise from different entities, and 
their behavior should not correlate, as is found for SN 2006tf. 
This model considerably diminishes the viability of the LBV 
hypothesis. 

In the above cases, and others in the literature, the high 
mass-loss rate was the primary factor in suggesting an LBV 
progenitor. We have shown that the high mass-loss rates 
may be incorrectly interpreted. This does not mean that 
the density is not high - a high density can be achieved 
independently of a high mass-loss rate. In fact our main 
point is that a high density should not be used as a proxy 
for a high mass-loss rate, and should not therefore be seen 
as indicative of an LBV progenitor. This is elaborated on 
further in fj3] 

2.4 Circumstellar Nebulae 

The similarity in the ring structures of the nebula around 
the LBV star HD 168625, the B supergiant Sh er 25, the star 
SBW1, and the nebula around SN 1987A led ISmithl (|2007f ) 
to argue that the progenitor of SN 1987A may also have been 
an LBV. While the first object is known to be an LBV, the 
next two have abundances inconsistent with a RSG phase 
ljSmithll2008bh 

An argument based purely on the shape of the CS 
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structures is quite misleading. It is true that the shapes 
of nebulae around LBV stars and around SN 1987A are 
quite similar. However, it should be noted that the these 
shapes are also similar to that of several bipolar plan- 
etary nebulae, which form around low mass starqj (< 
8Mq). Just as SNe cannot be attributed to be arising 
from low mass progenitors on the basis of their circum- 
stellar nebulae, they cannot be attributed to LBV progen- 
itors purely on the basis of the shaping of their CS struc- 
tures. The similarity in the shapes is a result of the sim- 
ilarity in the shaping mechanism, which is generally at- 
tributed to the interacting winds model and its variations 
jKwok et al.ll 19781; iFrank fc Mellemalll994l;lbwarkadas et all 
Il996l ; IDwarkadas fc Balickl Il998l ; iBalick fc Frank! |2002l l. "or 



2.5 SN 2005gj 



due to an inherent asymmetry in the wind s themselves 
|Frank et al.ll 19981 ; IDwarkadas fc Owockill2002l 'l. The nebula 
around SN 1987A has been modelled in a manner similar to 
the models for PNe ijBlondin fc Lundavisdll993h . The de- 
tails of the shaping mechanism may depend on the param- 
eters, but the basic ideas are similar. This does not imply 
that the stars themselves are similar. Thus the similarity in 
shapes of the surrounding medium cannot be argued as a 
similarity in the progenitor star. 

It could perhaps be argued that while the shaping mech- 
anism may be effective, in the interacting winds model it 
does take a change in wind parameters to bring about this 
effect. In the case of SN 1987A, this requirement was ful- 
filled by the star evolving from the blue to the red side of 
the HR diagram and then back to the blue before explod- 
ing. The change in wind parameters, indeed the entire ring 
system, has been at tributed to the transition from the RSG 
to the BSG stage l|McCravl [2007.. and references within). 
If such red-blue loops were always required for the forma- 
tion of the ring-like structures, this would set a timescale 
requirement which may be difficult to fulfill. However, we 
point out that such loops are but one way of achieving 
a change in the wind parameters. It has been reported 
that RSG and hypergiant stars may also sh ow variabil- 
ity, r e lated to changes in the wind parameters (|Humphrevsl 
2008; Lcvesque 200 9]) , and may even undergo pulsations 
jYoon fc Cantielldl2O10h . as described below. IRC+10420 
and VY CM a both show comp lex and extensive circumstel- 
lar nebulae ([Humphreys! [2008). As discussed below, some 
stars may experience a brief Wolf-Rayet phase before ex- 
plosion. Stars in binary systems can experience mass-loss 
episodes due to the pull of the secon dary. Mass ive stars are 
known to have weak magnetic fields (|Puls| [2008). which can 
influence the wind, leading to disk formation and strong 
shocks, or at the very least density enhancemen ts towards 
the magnetic equator l|ud-Doula fc Owockil l2002'l. A spheri- 
cal wind blowing into such a medium may lead to SN 1987A- 
like ri ng structures, as modelled by iBlondin fc LundqvislJ 
(1993). And finally, a rotating star may create such shapes 
due to the latitudinal asymmetry of the wind from the star 
itself l|Dwarkadas fc Owockil [2002). Given our still incom- 
plete knowledge of stellar mass loss and evolution, it would 
be therefore misleading to attribute such shapes purely to 
LBVs without concrete evidence. 



2 see catalog and images at 

http://www.astro.washington.edu/users/balick/PNIC/ 



iTrundle et al.l (|2008T ) find multiple absorption component P- 
Cygni profiles of H and He in the spect rum of 2005g j. Thes e 
profiles have earlier been noted by iPrieto et alj l|2007h . 
ITrundle et al.l (|2008h interpreted these as indicative of the 
progenitor's mass-loss history. Since such profiles have pre- 
viously been seen only in LBVs, the natural assumption is 
to connect SN 2005gj to an LBV progenitor. The similarities 
in the wind velocity deduced from the profiles and those of 
LBVs was used to further support their argument. 

The main problem with an LBV progenitor for SN 
2005gj is, that in the papers describing its disc overy and 
evolution (|Aldering et al.ll2006l ; [Prieto et al.ll2007i'). the type 
attributed to the SN was that of Type Ia. lPrieto et all (|2007f) 
did a cross-comparison with the SNID database, and the 
best comparison was overwhelmingly with a Type la spec- 
trum. It is clear that SN 2005gj is not a typical Type la, 
since it shows Balmer li ne emission, whose w idth seems to 
be increasing with time. ITrundle et alj (|2008f ) appropriately 
note that the typical S II and Si II, which solidified the 
la status of the very similar SN 2002ic (since classified as a 
Ia/IIn), are not as clearly delineated in SN 2005gj. However, 
the overall resemblance of the spectrum of SN 2005gj to that 
of SN 2002ic , betw een days 26 and 84, is clearly shown in 
IPrieto et all (120071'). The la-like featu res of SN 2002ic are 
even accepted by ITrundle et ail l|2008l ). Given this close re- 
semblance between the two SNe, it seems uncharacteristic 
that one would be a la and the other have an LBV progen- 
itor. If it is a core-collapse SN, it may offer an explanation 
for the absorption components, but the spectrum becomes 
more difficult to explain. 

The argument therefore boils down to whether SN 
2005gj is a Type la SN or not. If it is then it cannot have 
a massive star LBV progenitor under any circumstances, 
and that assumption is not tenable. It is possible that the 
multiple absorption components rightly signify mass-loss, 
but possibly from a companion star, not the progenitor. 
Type la's are presumed to occur in a binary system, so 
the fact that material from the companion star may sur- 
round the progenitor is not unexpected. This could, and has 
been used to, explai n the Balmer line emission in a Type la 
jHamuv et al.ll2005t . 

If it is shown to not be of Type la, then the LBV theory 
becomes more plausible. But it still needs direct proof that 
the profiles are due to the progenitor star as opposed to a 
companion. Also, it has to be shown that an LBV was actu- 
ally the progenitor, as opposed to the fact that the progen- 
itor star perhaps went through an LBV phase before losing 
its H and He envelope and undergoing a brief Wolf-Rayet 
phase (see jpp bef ore explosion. The spectra presented by 
IPrieto et aiTifeoOTh do show some resemblance to a Type Ic 
spectrum at late times. Furthermore, one then has to ex- 
plain the resemblance of the spectrum to that of other type 
la SNe, or at least to SN 2002ic, which would be difficult for 
a core-collapse event. Therefore, in the absence of further 
information, we remain agnostic towards the LBV interpre- 
tation. 



8 Dwarkadas 



3 DISCUSSION 

We have discussed in this paper various arguments put for- 
ward to suggest that luminous blue variable stars are pro- 
genitors of some core-collapse SNe. In many cases the high 
mass-loss rate that led to the suggestion of an LBV pro- 
genitor was based on the assumption that the medium into 
which the SN shock was expanding had a density profile de- 
creasing as r~ 2 , but the results show that the wind mass-loss 
parameter is variable with time, implying that the density 
profile does not decrease as r -2 but is more complicated. Ap- 
plying an equation that works for only constant mass-loss 
rate and velocity winds to every possible situation, without 
evaluating whether the inherent assumptions are satisfied, 
leads to incorrect results. Furthermore, unless the density 
profile is decreasing even more steeply than r -2 , application 
of these equations will generally yield a mass-loss rate that 
is higher than the true mass-loss rate. 

The case of SN 2005gl is an exception. The lack of a 
broad emission line in the spectrum 8 days after discovery 
distinguishes it from many of the other Type Iln. We have 
examined quantitatively several aspects of the arguments 
made by GL09, confirming the basic structure of a high- 
density medium surrounding the star. The identification of 
the progenitor star as having an exceptionally high luminos- 
ity in the V band (or close to it) is another argument in favor 
of an LBV candidate. However many questions still remain, 
such as whether an LBV-like star is the sole candidate for 
such a high-luminosity object, and why an LBV phase could 
have a high but steady mass-loss rate, resemblant of an LBV 
during eruption, for about 6 months, and then immediately 
explode as a SN. Given the currently available information, 
this could be considered a viable case for an LBV progenitor, 
with some unanswered questions remaining 

The case for quasi-periodic modulations in the ra- 
dio data is based on attributing those variations to 
S Doradus type instabilities. The arguments presented 
in this paper, when c ombin ed with those expressed in 
IChevalier fc Soderberd l|201(t ). seem to contradict the sug- 
gestion of an LBV progenitor here from several different 
points of view. However more detailed models and further 
observations are required for a conclusive assesment. The 
morphological similarity of nebulae around LBV stars to 
those around SNe such as SN 1 987A in no way implies sim- 
ilarity of the stars themselves, merely that of the shaping 
mechanism. 

The case of SN 2005gj presents a quandry. It could be 
considered a plausible case for an LBV, with reservations, if 
it can be conclusively shown that the SN was not of Type 
la, but is not tenable if the Type la interpretation is true. 
The jury is still out on this one. 

In summary, we have studied several suggested cases 
for LBV SN progenitors, and found that in only one case is 
there significant evidence to satisfy the interpretation of an 
LBV-like progenitor. While the difficulty in doing so may 
be considerable, it should be noted that in no case has vari- 
ability, a hallmark of an LBV, actually been demonstrated, 
although it has been alluded to. We emphasize here that 
this is not an argument for or against the fact that stars can 
explode in the LBV stage, and thereby act as progenitors of 
SNe. That is a discussion beyond the scope of this paper, 
and best left to stellar evolution theorists who are better 



equipped to answer that particular question. What we find 
here is that in only one proposed scenario is a strong case 
made to invoke an LB V-like object as the immediate progen- 
itor. 

Given our incomplete knowledge of stellar evolution, it 
is important to determine in the case of SN 2005gl whether 
other astronomical objects can have the required high lumi- 
nosity in the appropriate band, or whether an LBV is the 
only candidate. The other arguments, although consistent 
and supportive of an LBV interpretation, are not conclusive 
on their own. Even one LBV progenitor, if confirmed, does 
call for a modification of stellar evolution theory. The chal- 
lenge for stellar evolution theorists will then become more 
difficult, to explain not only why stars explode in the LBV 
stage, but whether they can do so immediately after a high- 
mass loss rate phase. 

At the same time, it appears that an LBV progenitor is 
not as pervasive as recent literature tends to suggest. It is 
clear that there are several Type Iln for which LBV progeni- 
tors, although suggested, are not necessary. This diversity in 
Iln progenitors provides some clues to their understanding, 
as discussed below. 

3.1 High Density CSM 

In many individual cases, the assertion of an LBV as a SN 
progenitor is made by noting that the observations require a 
high density CSM, and then assuming that the high density 
equates to a high mass-loss rate for the progenitor, typi- 
cally > lO _2 M0yr _1 , immediately preceding the SN explo- 
sion. The high mass-loss rate is attributed to an r\ Carinae 
type LBV explosion. It must be noted however, that even 
amongst LBVs, which are not very common in the first place, 
r) Car type explosive ejections are quite rare, with only two 
known events in the Galaxy - r\ car itself in the 1840s, and P 
Cygni. In these extreme events, the star actually increases 
in luminosity during the outburst. In other LBVs with S 
Doradus type variations, which comprise the more common 
variety, the bolom etric luminosity of the star does not vary 
(jPuls et al.ll2008T ). It remains to be seen if the frequency of 
explosive-type LBV events is high enough to account for all 
the SNe that are attributed to it. 

Although we have shown that in many cases the high 
mass-loss rate argument may not hold, the high X-ray, radio 
or optical flux suggests that in many scenarios the CSM den- 
sity close to the star must be high. It is therefore prudent to 
mention here that high density does not necessarily imply a 
high mass-loss rate. Other mechanisms can result in a high 
density medium. A common mechanism is the interaction of 
the fast wind from the star with the slower wind from a pre- 
vious epoch, forming a wind-blown bubble. The fast wind 
sweeps the ambient medium or slower wind into a dense 
shell, whose density can be extremely high if the wind shock 
is radiative (see Figure 1), and the swept-up mass collapses 
into a very thin dense shell bordere d by a radiative shock 
llWeaver et al ]| 19771 ; |Pwarkadasll2005l ) . The shell density and 
thickness depend on the exact wind parameters and inter- 
action time, but the density can exceed 10 4 cm -3 . The mass 
in the dense shell comes from the swept-up material, not 
the progenitor wind, and could have been emitted over tens 
of thousands of years. Depending on the ambient density, 
the wind parameters and the time taken for the interaction, 
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the size of the bubbles and the density enhancements may 
vary. A wind expanding supersonically for a sufficiently long 
time can lead to such dense shells, more-or-less independent 
of the mass-loss rate. The mass-loss rates cannot then be 
obtained from the density using the formulae described in 
1231 

It is well known that as stars evolve, their wind param- 
eters change, lea ding to the formation of large wind bubbles 
around the star (iGarcia-Segura et al.lll996l: ICTiuI|2003|, 120081 ; 
iFrever et al, I l200l [2006: Arthui] l2007l ; iDwa rkadas 2007alldh. 
Moderate changes in the wind parameters lead to the forma- 
tion of smaller versions of these bubbles during the last few 
thousand years of a star's life. SN 1987A for example is sur- 
rounded by a win d bubble forme d about 20,000 years before 
the SN explosion (|McCravH200"7h . The progenitor's mass-loss 
rate must have been quite low , not exceeding 10~ 8 Me ) yr~ 1 
JChevalier fc Dwarka das 1995 1; lLundqvistlll999l ; Dwarkadasl 



l2007d lbh. The calculations of [Dwarkadas et all l^Old ) infer 
the presence of a wind bubble around SN 1996cr, surrounded 
by a dense shell at about 0.03 pc. 

Figure [T] shows the structure of such a bubble around 
a massive star, taken from a simulation of a Wolf-Rayet 
wind with mass-loss rate of 2.5 xlO _6 MQyr _1 interacting 
with a RSG wind from a previous epoch (see for example 
I Dwark adas 2007d|). Note that the density contrast between 
the density inside of the bubble and the shell is more than 
five orders of magnitude. The wind termination shock (Rt), 
the contact discontinuity between the two winds (Rcd) and 
the outer shock (R a ) are marked. Significantly, in this struc- 
ture, only the region labelled as free wind (to the left of the 
wind termination shock Rt and the dotted line) has a den- 
sity profile that decreases as r -2 . Equations [1] to [3] above can 
only be applied to this region to compute the wind parame- 
ters. If these equations are used to compute the density at 
other radii, especially in the dense shell region, it will result 
in a derived wind parameter that is extremely high, in a case 
where there is no freely expanding wind. This can mislead 
one into believing that the progenitor is a high mass-loss rate 
LBV star. The situation in a star with a varying mass-loss 
rate and wind velocity will be even more complicated. 

A high density can also be achieved in localized regions 
due to high density clumps, which may be formed by hy- 
drodynamical instabilities in the pre-existing circumstellar 
medium or in the ejected material. The interaction of the SN 
blast wave with a clump leads to a slow-moving shock ex- 
panding into the clum p, and a reflected shock back into the 
inter-clump medium l|Klein et al.l 1 19941 ) . The slow- moving 
shock velocity is represented in the formation of narrow lines 
(see i]3.3[) . The clumps, although high density and sufficient 
to give rise to the narrow lines in the spectrum, make only 
a minimal contribution to the total mass of the medium, 
are not representative of the overall density of the medium, 
and cannot be used to derive the global properties of the 
medium. 

Finally, as current stellar evolution theor y suggests for 
stars that are more massive than about 40 Mq (jHirschi et al.l 
|2010T ). the possibility exists that the massive star experi- 
enced an LBV phase before becoming a Wolf-Rayet star. 
If the Wolf-Rayet stage is short enough, the dense LBV 
wind may be close in to the star, and/or may be swept- 
up by the faster W-R wind. Either of these possibilities 
will give rise to denser material in the surrounding medium, 



due to an LBV phase, but does not require the LBV to 
be the progenitor. This kind of s cenario was suggested 
fo r SN 2001em by ICheyaliei) l|2007h . following the models 
of IChugai fc Chevalierl (120061 ). which showed that a dense 
medium around the SN may have formed about 2000 years 
earlier, with a mass-loss rate on order 10 -3 MQyr -1 . The star 
may have transitioned from an LBV to a W-R star shortly 
before explosion. 



3.2 Short Wolf-Rayet Phase 

In the case of SN 200 lem , SN 1996cr, and many oth- 
ers ((Dwarkada s et al]|2010l . and references therein; also see 
i]3.4p . models suggest that the star explodes after a short 
W-R phase (< 10 4 years). Stellar evolution models, on the 
other hand, generally show that the W-R phase lasts for > 
10 years. Under what conditions could we reasonably ex- 
pect the star to remain in the W-R phase for only a few 
thousand years or less before explosion? It turns out that 
there are several possibilities. (1) One way this could hap- 
pen is if the star has a binary companion. lEldridge et al.l 
(2008) have shown that stars in a binary system may en- 
ter the W-R phase at much lower initial mass as compared 
to single stars. The stars at the very low end of the mass 
range of stars able to become W-R sta rs will have very shor t 
W-R lifetimes (see Figures 2 and 3 in lEldridge et"aill2008l '). 
on order tens of thousands of years or less. In principle, the 
same short lifetime would apply to single stars which are 
at the lower end of the mass range where they can become 
W-R stars, although since the cutoff mass to become a W-R 
is higher, the number of stars is much smaller. (2) Binary 
interaction can result in Case C mass transfer which follows 
the He burning phase. Under certain cond itions this can 
lead to explosive common-enve lope ejection (|Podsiadlowskil 
120071 ; iPodsiadlowski et alj|2oioh , which leads to the ejection 
of both the hydrogen and helium layers late in the evolu- 
tion of the star. This results in a W-R star with a natural 
10 4 year or less timescale. (3) Another promising channel for 
the formation of a sh ort duration W-R s tage w as recently 
shown by the work of lYoon fc Cantieliol J20 ldf) , which fol- 
lowed up on earlier work by Heger et al.l (|l997T ). They find 
that RSGs could have strong pulsation driven superwinds. 
These would cause the RSG to just remove its outer layers 
and become a W-R star in the last hundreds to thousands of 
years of evolution. Such pulsation-driven winds may also re- 
sult in the formation of circumstellar shells around the star, 
and thus may be usef ul in understanding Type Ilns. In fact, 
iFransson etaH (|2002T ) suggested that Type Iln SNe, and es- 
pecially SN 1995N, could arise from RSG progenitors which 
experience superwinds, a scenario that resembles the one 
outlined above (although they did not specifically advocate 
a W-R progenitor). 

Overall it seems that several channels exist which can 
result in short-duration Wolf-Rayet phases, and therefore 
it is reasonable to propose that some SNe can experience 
a short W-R phase before explosion. Many scenarios that 
require a short W-R phase tend to indicate stars with ini- 
tial mass that is lower than the minimum single-star mass 
required to become a W-R star. If some fraction of Type 
Iln arise from such lower mass W-R stars, then it may in- 
dicate that, contrary to the general assumption, Ilns do not 
arise from very massive stars but instead from a lower av- 
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Figure 1. The density profile in a wind-blown bubble around a high-mass star. The bubble was formed by a Wolf-Rayet wind, with 
a mass-loss rate of 2.5 xl0 — 6 M0yr _1 and a wind velocity of 2500 km s — 1 , interacting with a RSG wind from a previous epoch. The 
interaction gives rise to a low-density bubble surrounded by a high density shell. This is one way to get a high density in the surrounding 
medium without have an extremely high mass-loss rate wind. In this figure, the region marked as "Free Wind" denotes the freely 
expanding wind, to which equations [T] to [3] can be applied. Application of these equations at other radii could result in very large, but 
incorrect, mass-loss rates. 



er a^e initial mass pop u lation . A similar result was derived 
by Anderson fc James] ((2008), who found, from an investi- 
gation of the association of the explosion sites of Type Iln 
SNe with recent star formation, that the majority of Type 
Ilns do not arise from the most massive stars, which would 
be the case if they all had LBV progenitors. 



Blue supergiants (BSGs) are also known to have velocities 
of a few hundred km s _1 . The circumstellar medium close 
in to SN 1987A has been mod elled with a BSG wind ve- 
locity between 300-600 km s~ 1 (iBlondin fc Lundavist]|l993l ; 
IChevalier fc Dwarkadail 19951 ; |Pwarkadasll2007d l 



3.3 Velocities indicated by Narrow Emission lines 

In many of the SNe discussed, and several others, nar- 
row lines seen in the spectrum are invoked as arising from 
the unshocked, ionized CSM. In some cases, the presence 
of P-Cygni absorption components may indicate the pres- 
ence of an outflow. Their velocity is interpreted as the 
velocity of the ambient wind, and is considered interme- 
diate between that of a RSG an d a W-R but appropri- 
ate for LBVs (|Trundle et al.ll2008h . However , it is not cer- 
tain w hat appropriate velocities are for LBVs. [Trundle et al.l 
il2008h state that LB V velocities lie between 50-300 km s" 1 , 
iKotak fc Vinkl d2006h sugge st velocities between 100 and 500 
km s _1 , (|Smith et alJbOlCO find 600 km s _1 from their spec- 
trum of SN 2005gi appropriate for an LBV, and r \ Cari- 
nae ha s reported velocities exceeding 3000 km s _1 l|Smithl 
2008a). Furthermore, such velocities are not unique to LBVs. 



As discussed above there are other equally plausible al- 
ternati ves for the na rrow emissi on lines. In the cases of SN 
1986J dChugailll993h, SN 1988 z (|Chugai fc Danziger|[l99l l. 
SN 1978K (jchugai et al.|[l995h and SN 1996cr (Dwarkadas 
et al. 2010) at least, lines of narrow or intermediate width 
are explained as arising from the interaction of the SN shock 
with a clumpy medium. The FWHM of these lines denotes 
the velocity of the shock driven into the clumps. They are 
unrelated to the velocity of the surrounding medium. There 
also exists the possibility that the narrow line spectrum is 
contaminated by emission from the host galaxy. We note 
that clumps could also provide the high density for much of 
the emission. Interestingly, in the case of SN 1996cr, both 
a dense shell (to provide the enhanced emission) and dense 
clumps in the ambient wind (to explain the na rrow lines) 
have been postulated bv lDwarkadas et al.l (|2010t ). 
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3.4 Type Iln SNe 

The Type Iln SN class is most likely to be associated with 
LBV progenitors. However Type Iln's are not a homoge- 
neous gro up, and in fact not all of them show strong radio 
emission l|van Dvk et al.l fl996) or X-ray emission. This may 
suggest that not all Type Iln are undergoing CS interaction 
with a high density medium. The diversity of Ilns makes 
it inappropriate to classify Ilns, or the similarly designated 
Type Ibn or Icn, as a separate class of SNe. Some are more 
indicative of a phase in the lifetime of a SN, where it is inter- 
acting with high density material, and which it will outgrow 
in a matter of years or decades, no matter what class they 
fall in or the nature of their progenitor. A similar point had 
been made earlier bv lKotak et al.l l|2004 ) with regards to the 
dual classification of SN 2002ic as a Ia/IIn. The fact that 
Type la's sometimes share a dual classification as Type Ilns 
typifies the inhomogeneity among the group, and their lack 
of a common progenitor. 

Many Ilns earn the 'n' designation later on in their life- 
time, when they are presumably interacting with high den- 
sity material and are visible. But this merely emphasizes 
a feature in the late-time spectrum rather than a physical 
property of the progenitor. Mainly it is the initial designa- 
tion of type which would be useful in pinpointing a progen- 
itor. This is a small but subtle distinction, because unlike 
other SNe classes, it suggests that Type Iln's need not arise 
from a single progenitor (if in fact other subtypes do arise 
from a single progenitor, which has not been shown). The 
fact that Ilns could have both high and low initial mass 
progenitors is consistent with this statement. 

A scenario of a wind bubble blown by a W-R or blue 
su pergiant wind was sugge sted for the Type Iln SN 1996cr 
bv lDwarkadas et al.l l|201Cf ). with the mass- loss rate of the 
progenitor star being < lO -4 M0yr -1 , and in all probability 
su bstantially le s s. Th e Type Ibn SN 2006jc was suggested 
bv iFolev et al] |Io07|) to be a Wolf-Rayet star exploding 
in a dense He-rich m edium formed about two years be- 
fore the star exploded. iPastoreilo et al l (|2008l ) propose that 
the progenitor of the hybrid SN 20051a was a very young 
Wolf-Rayet (WN-type) star which experienced mass ejection 
episodes shortly before core collapse. SN 2003bg was initially 
classified as a Ic before becoming a Iln, while SN 2001em 
transformed from a Ib/c to a Iln, possibly undergoing sub- 
stantial mass-loss about 1000-2000 years before explosion 
jChugai fc Chevalieij[200ej ). IStockdale et all l|20ld ) indicate 
that SNe 1996aq and 2004dk, both classified as Type Ic af- 
ter explosion, are evolving in a similar manner to SN 1996cr 
and SN 2001em. 

These considerations suggest that a subset of SNe with 
the 'n' designation may show a wind-bubble structure on 
parsec or even sub-parsec scales, due to a change in the 
wind parameters close to the end of the star's life. We pro- 
pose that one channel for the formation of Type Ilns that 
show enhanced late-time X-ray and radio luminosity may be 
W-R stars with fast winds, where the final W-R phase lasts 
for a few thousand years or less (see t)3.2p . The W-R wind 
sweeps up the surrounding medium, leading to the formation 
of dense shells close in to the star. The interaction of the SN 
shock wave with this medium gives rise to many of the no- 
table properties of Type Iln SNe, including the large X-ray 
and radio luminosity at late times, while obviating the need 



for a high mass-loss rate progenitor (the mass-loss rates f or 
Galactic W-R stars are < lO-^Moyr" 1 l|Crowtheijl26og )1. 

There are some Type Ilns that show Iln features al- 
most from the ti me of explosion, including SN 2005i p 
jSmith et all 120091 1 and SN 1995N jFransson et all 120021 ). 
In these SNe it appears that a high density region sits just 
outside the stellar envelope. SN 2005gl also shows an ex- 
tremely high density region outside the stellar envelope, but 
as seen it does not extend far out radially. These are all 
indicative of a wide diversity in this class of objects, and 
would probably require a different channel for their forma- 
tion. SN 2005gl has been associated with an LBV-lik e pro- 
genitor. Using Spitzer spectroscopy, I Fox et all (|2010l ) have 
studied SN 2005ip, and conclude that a progenitor erup- 
tion formed a dust shell about 100 years prior to explosion. 
They find that a large mass-loss rate > 10~ 2 M© yr^ 1 is re- 
quired to explain this dust shell, and therefore suggest that 
it too is consistent with an LBV progenitor. On the other 
hand, ISmith et al l {2009) found a mass-loss rate 2-3 orders 
of magnitude lower. Their optical observations do not re- 
quire, or invoke, a dust shell, and they assumed that the 
mass lo ss arises in a stead y wind with density decreasing 
as r " 2 . ISmith et all l|2009h concluded that the progenitor 
was a RSG star, despite having a wind velocity of about 
120 km s _1 . It is interesting to speculate whether these two 
arguments can perhaps be r econciled by invoking pu lsation- 
driven superwinds in RSGs (|Yoon fc Ca nticllo 201CT), which 
could potentially give rise to the high mass-loss rate dust 
shell just prior to explosion. 

Follow-up multi-wavelength observations with increased 
frequency, specifically in the radio and X-ray regimes, will 
undoubtedly reveal more Type Iln SNe. Regular monitor- 
ing may detect more analogs of SN 1996cr, with X-ray and 
radio emission increasing over year-long periods of time. 
Such observations, accompanied by detailed theoretical in- 
vestigations including hydrodynamic modelling and multi- 
wavelength emission computations, are urgently needed, if 
we are to determine the nature and progenitors of this in- 
teresting category of objects. 
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